About possible contribution of intrinsic charm 
component to inclusive spectra of charmed mesons 



E. V. Bugaev and P. A. Klimai 

Institute for Nuclear Research, Russian Academy of Sciences, 60th October 
Anniversary Prospect 7a, 117312 Moscow, Russia 

E-mail: pklimai@gmail.com 

Abstract. We calculate differential energy spectra (a^ -distributions) of charmed 
particles produced in proton-nucleus collisions, assuming the existence of intrinsic 
heavy quark components in the proton wave function. For the calculation, the recently 
proposed factorization scheme is used, based on the Color Glass Condensate theory 
and specially suited for predictions of a production of particles with large rapidities. 
It is argued that the intrinsic charm component can, if it exists, dominate in a sum of 
two components, intrinsic + extrinsic, of the inclusive spectrum of charmed particles 
produced in proton- nucleus collisions at high energies, in the region of medium xf, 
0.15 < xp < 0.7, and can give noticeable contribution to atmospheric fluxes of prompt 
muons and neutrinos. 
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1. Introduction 

According to the usual assumption, heavy quark partons in the nucleon are radiatively 
generated, i.e., they appear only as a result of QCD evolution which starts from a 
null value at a scale of approximately the corresponding quark mass. There are two 
reasons, at least, for this assumption: i) if quark mass is large enough, the partonic 
probability distribution function (PDF) is perturbatively calculable and ii) up to now, 
exact experimental constraints on these PDFs are absent. 

The mass of charm quark (~ 1.3 GeV) lies in between the soft and hard scales, 
so, the standard method of PDF's calculation is not well justified. Besides, from a 
theoretical point of view the existence of an "intrinsic charm" (IC) component of the 
nucleon is not forbidden. Many nonperturbative models (the light-cone models [TJ [2] , 
the meson cloud picture j3j Hj, the MIT bag model [5], SU(4) quark model [6], etc.) 
predict such a component at an energy scale comparable to mass of charm quark m c . 

In the present paper we use the predictions of the specific light-cone model 
developed in [TJ ("BHPS model"). In [1] the PDF of charm quarks from the 5-quark 
component uudcc in the proton is calculated. The uudcc- Fock state arises in QCD not 
only from gluon splitting which is included in DGLAP evolution, but also from diagrams 
in which the heavy quark pair is multi-connected to the valence constituents. The latter 
components are called intrinsic charm Fock components. It was originally suggested in 
[1] that there is ~ 1% probability of IC Fock states in the nucleon. The materialization 
of the IC Fock state (after its interaction with the target) leads to the production of 
open-charm states such as scalar mesons D + (cd), D°(cu), or vector mesons D* + at large 
values of Feynman variable xp- This may happen either through the coalescence of 
the valence and heavy quarks of the same rapidity, or through hadronization of the 
produced heavy quarks. In the present paper we do not consider the process of the 
coalescence (and the corresponding leading particle effect) and assume that quarks of 
cc-pair fragment into hadrons independently of each other. 

The x-dependence for the PDF predicted by the BHPS model is 



Here, /io is the starting value for the factorization scale, /io = m c . The normalization 
constant A is a variable parameter of the model. The magnitude of IC can be 
characterized by the momentum fraction carried by c + c, 



The global QCD analysis [9] rules out the possibility of an IC component with the 
momentum fraction (x) c+5 much larger 0.02 (for fi = fi = m c ). 

From a point of view of QCD, it is natural to expect that a contribution of heavy 
quark component in hadron's wave function cannot be large, due to a suppression by 
the power of heavy quark mass squared. However, it is possible that this suppression is 
compensated by sufficiently strong nonperturbative gluon fields in the nucleon (see, e.g., 



f c (x,fi ) = f £ (x,fi ) = Ax 2 6x(l + a;)lna; + (1 - x)(l + lOx + x 2 ) . (1) 



i 




(2) 



o 
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[7]). It had been shown in [5] using the technique of heavy quark mass expansion, that 
the intrinsic charm contribution to the nucleon momentum, (x) c+5 , at the scale \i ~ m c , 
is proportional to the (unknown) factor A (which is determined by the matrix elements 
of gluon operators) whose value can be obtained using nonperturbative methods in QCD 
(e.g., a theory of instanton vacuum). If A is of order of typical strong interaction scale, 
A ~ 1 GeV, then, according to [Sj, the prediction is 

(x) c+ - c = few x 1(T 3 , (3) 

that is quite similar with the suggestion of pp. 

We consider also, as a second example, the phenomenological model [9], in which 
the shape of the charm PDF is sea-like, i.e., similar to that of the light flavor sea quarks 
(except for normalization). Namely, for this case [9] 

f c (x, no) = f s (x, no) ~ fe(x, m c ) + f a (x, m c ). (4) 

It is very important, for high energy cosmic ray physics, that the heavy constituents 
of an intrinsic heavy quark state, such as \uudcc), usually carry the largest fraction of 
the nucleon's momentum. This can easily be shown in the context of light-cone models, 
where contributions of high-mass Fock states (to probability distributions, for example) 
are suppressed only by energy denominator factor [TO] 

-2 

(5) 

Here, j is a constituent number, Xj = Pj + ^ /p^ are the light-cone momentum fractions, 
= (p(°> + p^)/\/2, p^ is a particle's energy and p^ is its longitudinal momentum 
(i.e., the momentum in a direction of the parent nucleon). Evidently, if rrij is large, the 
value of Xj should also be large, for far off-shell configurations not to be too suppressed. 

In two recent works [Hi [12] the problem of a search of intrinsic charm was again 
considered in connection with experiments at RHIC, LHC and Tevatron. These authors 
argue that the hypothesis of intrinsic charm can be checked in these experiments 
(studying the resulting transverse momentum spectra of charmed mesons). In the 
present paper we calculate inclusive distributions of charmed mesons (a^-spectra) 
produced in p — collisions at high energies, using the mechanism proposed in [TT] 
(based on the CGC formalism). 

2. Hadron production in the forward region and CGC 

We are interested in heavy meson production in p/l-collisions at very high energies and, 
what is most essential, at the projectile fragmentation region, i.e., at forward rapidities. 
In the forward rapidity region the Bjorken x of the target nucleus has its lowest possible 
value for a given a/s whereas the Bjorken x of the projectile is large (close to unity). 
Rescatterings of the projectile's quarks in a target are quite substantial. It had been 
shown in [13], [H] that a particle production in p-nucleus collisions at high energies and 
forward rapidities can be expressed through the cross section of a multiple scattering 



S 



" p 2 j± + m) 

Z-~t rr- . 

3=1 X 3 
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of a quark from the target nucleus which is treated as a Color Glass Condensate. In 
this approach the projectile proton is considered as a collection of quarks and gluons 
(according to the parton model), and the standard DGLAP evolution and collinear 
factorization are employed. The resulting formula [15] for a single- inclusive charm 
production in pA collision is 

r/A (pA -> DX) = -J— [ dx p ^ x 



dyd 2 p± (2tt) 2 J xp 



Xp 



X r 



xfc/ P (x P ,Q 2 )N F (^p ± ,x A )D D/c (^,Q 2 ). (6) 

Here, p±, y, xf are the transverse momentum, rapidity and Feynman x of the produced 
D meson, respectively. The variable x p is a momentum fraction of the charm quark of 
the projectile, xa is a momentum fraction of the target parton. The connection between 
xa and x p is 

x A = ^M±Z4j = xe -2v (7) 
xps 

and the rapidity is connected with xp, p± by the formula 



x F = = e y . (8) 

One can see from ([7]) and (jSJ) that at large rapidities the x^-values are small while x p is 
large. 

The functions f c /p{x p ,Q 2 ) and Dp>/ c (z,Q 2 ) are the proton charm PDF and the 
charm fragmentation function into D mesons, respectively. The factorization scale Q 2 
is taken to be equal to p\. At last, the function Np(q±,XA) is equal to an imaginary 
part of the quark- ant iquark dipole- nucleus forward scattering amplitude, q± = ^P±- 
This amplitude can be modelled by the following expression [16J: 

1 f •? ^9/ ,\~y(q±,x A ) 



N F ( q±J x A ) = J d 2 r ± e^- [l - exp (-i (t^xa)) 1 ^) 



(9) 



Ht) Xa "> (10) 



Here, r± is the transverse size of the dipole, Q s (x) is the saturation scale (see, e.g., [T7] 
for a definition of Q s ), parameterized by the simple formula 

Ql(x) N 2 - 1 
1 GeV 2 2N 2 

x = 3 x 10~ 4 , A = 0.3, N c is a number of colours (iV c = 3). For the "anomalous 
dimension" factor 7((?_l,:ea) (which enters Eq. Qj due to quantum evolution in the 
Color Glass Condensate), we use the parametrization from [IS] : 

w a - 1 

7M=-ft + (i-»0 (m ._ 1) + t . (ID 



where 
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and a = 2.82, b = 168, 71 = 0.628. These parameters were fitted to the data of the 
RHIC experiment [19] where deuteron-gold collisions at forward rapidity (i.e., at the 
fragmentation region of the projectile) were studied. One should note that just this 
region of rapidities is essential for calculations in cosmic ray physics. 

The formula (|6]) gives the spectrum of produced particles normalized on their 
average multiplicity in one collision. Correspondingly, on vertical axes of Figs. [T] - 
|3]we present numbers of produced particles in the p^4-collision, in unit interval of xp. 

For a calculation of the inclusive cross section one must multiply dN in (j6]) on a 
cross section of the collision, i.e., on a transverse size of the target nucleus. To determine 
an approximate value of the proportionality coefficient, we note, firstly, that the cross 
section formula corresponding to Eq. (jSJ) is given by (151 [20] 



da(pA -»■ DX) : ' 



/dxp fc/p(*Epi Q 
Xf 



^ X 



dyd 2 p± (2vr) 2 J xp' 

' d 2 r ± e^ a J^llA . D D/C (^, Q 2 ). ( I 3 ) 



Xp 



where <Tdi P is the dipole-nucleus cross section which, by optical theorem, is expressed 
through an imaginary part of the dipole-target forward scattering amplitude (in r±- 
representation): 

ffdip (rj. , x) = 2N F (r ± ,x). (14) 

The Fourier-Bessel transformation relates Np(r±,x) with the amplitude in q±- 
representation, introduced above, in Eq. (J9]): 

J d 2 r ± e l<r ^Np{r ± ,x) = 

= 2n J rxdr ± J (q±r±)N F (r ± ,x) = N F (q ± ,x). (15) 

For the estimate of the proportionality coefficient, we use, for simplicity, the 
quasiclassical approximation, in which anomalous dimension factor, •y(q±, x), is equal to 
1. The approximate expression for the dipole cross section can be written as [2Tj 

■K 2 r 2 ± a s (ii 2 )xG A (x } /j, 2 Y 



(?dip(r±,x) = a 



1 — exp 



(16) 



3<T 

where [i 2 = /Xq + C/r\, Ga(x, fi 2 ) = AG(x,fi 2 ) is the gluon density in the nucleus, //q, 
C and do are parameters. The expression ( II 6p agrees, in the limit of small r±, with the 
leading twist perturbative QCD result. Using, in the region fi 2 ~ fj,Q, the definition of 
the saturation scale Q s pj 



Q 2 s (x) = 2n 2 ^a s { l il)^xG A {x ll il), (17) 

the exponent in ffl6|) can be rewritten in the form of (Q (with 7 — > 1), if Co = 2ttR\. 
One can easily check that the Eq. (|T7|) is not in contradiction with the parametrization 
of Eq. ( TlOl) because one has, approximately, 

xG A {x) ~ Ax-\ R\ ~ A 2/3 . (18) 
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Finally, it follows from these considerations that the coefficient, connecting da and 
dN, is approximately equal to ^ = ttR\. It is known (21] that, for a dipole-nucleon 
cross section, the value of in Eq. ( jl~6l) is equal to ~ 20 mb. This value must be 
scaled according to the transverse area of the target nucleus. For the nucleus one has 
Ra ~ 1.3v4 1//3 fm, and for the nucleon one has Rn ~ 0.8 fm. 

For the numerical calculations in the present paper we used the value of A equal to 
14 having in mind that this value is close to the average atomic number of the nuclei in 
atmosphere. For A = 14 the proportionality coefficient is 

^ = n R A « 300mb. (19) 
2 



3. Results and discussions 



For the concrete calculations we used the charm PDFs presented in [9]. Authors of [9] 
used for the presentation of PDFs at different values of the factorization scale the results 
of CTEQ group [221 |9]. For BHPS model, it was assumed, in the present paper, that 
the charm content of the proton is on the maximal level, {x) c+5 = 0.02. The same is for 
the sea-like IC model: (x) c+5 = 0.024. 

For applications in cosmic ray physics (in particular, for calculations of atmospheric 
spectra of leptons (see, e.g., [231 121])) we need the dependence of dN/dxp on xp, so, we 
integrate (EJ) on p± and change the variables, using the relation 

dxF (on\ 
— — = x F . (20) 

dy 

The results of the calculations are shown in figures [I] - HI 

The z distributions of the c-quark fragmentation functions D(z, Q 2 ) at their starting 
scales were assumed to obey the Bowler parametrization [25J, 

D(z,Ho) = N c z- (1+ ^\l - z) ac e-^ /z , (21) 

with three parameters N c , a c , 7 C . In the present paper we did not calculate the DGLAP 
evolution of fragmentation functions, and use f l2T|) for all values of /i. 

The concrete values of the parameters N c , a c , 7 C were taken from the work [25] , 
in which the nonperturbative fragmentation functions for D-mesons had been obtained 
by fitting experimental data on single-hadron inclusive production in e + e~-annihilation 
obtained by Belle, CLEO, ALEPH and OPAL collaborations. 

In figure [1] the inclusive spectra of D°-mesons for several values of cms-energy 
are shown (for BHPS model). It is seen that the steepness of the spectra strongly 
varies with xp: the relatively flat region ends at xp ~ 0.5 -j- 0.6. Figure [2j where the 
comparison of the BHPS and sea-like models results is given, shows that the form of xp- 
spectrum is strongly model-dependent. Figure |3] shows that the x^-spectra for different 
charmed mesons are quite similar: the slight difference is connected, evidently, with the 
corresponding difference in PDFs of light sea antiquarks in the nucleon. 

Inclusive cross sections of charmed mesons had been calculated in many works. We 
compare our results with the corresponding results of the most recent works [271 EH]- 
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Authors of [271 EH] calculated the cross sections of inclusive cc-pair production using 
the colour dipole-model approach [291 [301 [31] ( an d extrinsic charm). The results of [27] 
and |28j are quite similar; we compare our curves for E p = 10 9 GeV and E p = 3.5 x 10 4 
GeV (for BHPS model) with the corresponding curves of [27] (see figured]). For this 
comparison we multiply our curves for dN/dxF on 2 ■ ttR\ ~ 600 mb. The coefficient 2 
takes into account, approximately, the contribution of other D-mesons (D + and D* + ). 
Figure H] shows that the contribution of the intrinsic charm component can be relatively 
large and, in particular, sufficient for applications in cosmic-ray experiments. 

It follows from figure H] that in the tail region, xp > 0.7, our curves are even steeper 
than those of [27], in spite of their "intrinsic" origin. One must bear in mind, however, 
that in the present paper inclusive spectra of charmed hadrons rather than cc pairs 
are calculated. Naturally, a taking into account of the fragmentation of c quarks into 
hadrons leads to additional steepening of x^-spectra. 
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Figure 2. The comparison of Xi?-spectra for BHPS (solid lines) and sea-like (dashed 
lines) models, for two energies, y/s = 10 2 and 10 4 GeV (A = 14). 
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Figure 3. The comparison of rr^-spectra for different D-mesons, for y/s = 10 4 GeV 
and A = 14. Thick solid curve is for D°, dashed curve - for D + , thin solid curve is for 
D*+. 
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Figure 4. The comparison of the present results (solid lines) with the results of [27] 
(dashed lines), for E p = 10 9 GeV (upper curves) and 3.5 x 10 4 GeV (lower curves) 
(A =14). 



